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ABSTRACT 

Mapping Mg II resonance emission scattered by galactic winds offers a means to determine the spatial 
extent and density of the warm outflow. Using Keck/LRIS spectroscopy, we have resolved scattered 
Mg II emission to the east of 32016857, a star-forming galaxy at z = 0.9392 with an outflow. The 
Mg II emission from this galaxy exhibits a P-Cygni profile, extends further than both the continuum 
and [O II] emission along the eastern side of the slit, and has a constant Doppler shift along the 
slit which does not follow the velocity gradient of the nebular [O II] emission. Using the Sobolev 
approximation, we derive the density of Mg^ ions at a radius of 12 - 18 kpc in the outflow. We model 
the ionization correction and find that much of the outflowing Mg is in Mg^^. We estimate that the 
total mass flux could be as large as 330 - 500 M© yr~^, with the largest uncertainties coming from 
the depletion of Mg onto grains and the dumpiness of the warm outflow. We show that confining 
the warm clouds with a hot wind reduces the estimated mass flux of the warm outflow and indicates 
amass-loading factor near unity in the warm phase alone. Based on the high blue luminosities that 
distinguish 32016857 and TKRS 4389, described by Rubin et al. (2011), from other galaxies with 
P-Cygni emission, we suggest that, as sensitivity to diffuse emission improves, scattering halos may 
prove to be a generic property of star- forming galaxies at intermediate redshifts. 

Subject headings: 



1. INTRODUCTION 

In the baryon-driven picture of galaxy evolution, gas 
accretion onto galaxies, gas consumption by star forma- 
tion, and gas ejection by galactic winds shape the prop- 
erties of galaxies. The only simple aspect of this baryon 
cycle may be the rate of gas infall into galactic halos, 
which may follow the mean growth rate of the dark mat- 
ter. How much of the inflow can cool and fuel star forma- 
tion, versus being shock heated to the halo virial temper- 
ature, remains a largely theoretical debate (Keres et al. 
2005; Dekel & Birnboim 2006; Keres et al. 2009; Dekel 
et al. 2009; Bouche et al. 2010; van de Voort et al. 2011; 
Dave et al. 2012). In practice, it has generally been more 
straightforward to detect outflowing gas than inflowing 
gas. 

Over a broad redshift range, spectroscopic surveys have 
empirically described the galaxy population hosting out- 
flows using blueshifted, resonance absorption lines, where 
the sign of the Doppler shift (relative to the galactic 
barycenter) unambiguously identifies outflowing gas on 
the near side of a galaxy. (Shapley et al. 2003; Mar- 
tin 2005; Tremonti et al. 2007; Sato et al. 2009; Weiner 
et al. 2009; Chen et al. 2010; Steidel et al. 2010; Heck- 
man et al. 2011; Coil et al. 2011; Erb et al. 2012; Martin 
et al. 2012, Rubin etal. , in prep). These data provide 
an accurate picture of which galaxies host outflows, e.g., 
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namely those with large concentrations of massive stars 
(Heckman 2003; Kornei et al. 2012; Law et al. 2012) con- 
sistent with theoretical arguments requiring massive star 
clusters to generate outflows (Murray et al. 2011). 

The properties of these outflows, in contrast to their 
demographics, remain poorly constrained. The outflow- 
ing mass flux, for example, is central to galaxy evolu- 
tion models; yet theoretical models that fit the empiri- 
cally determined mass function and mass - metallicity 
relation for galaxies require a mass-loss rate of order 
the star formation rate (SFR) at mass scales compa- 
rable to the Milky Way (Dave et al. 2012; Shen et al. 
2012; Creasey et al. 2013; Hopkins et al. 2013; Puch- 
wein & Springel 2013). The location, maximum veloc- 
ity, and total column density of the warm outflowing gas 
are not well-constrained by observations of optical and 
ultraviolet resonance lines in galaxy spectra. The ab- 
sorbing clouds may lie anywhere along the sightline to 
the galaxy, and the decline in gas covering fraction with 
increasing blueshift makes it challenging to detect the 
highest velocity gas (Martin & Bouche 2009). The hot- 
ter wind fluid, which entrains this warm, low-ionization 
gas, largely eludes direct detection at intermediate red- 
shifts (Tripp et al. 2011). 

Two complementary strategies have emerged recently 
for measuring the spatial extent of low-ionization out- 
flowing gas. First, at large radii, outflows can be mapped 
with sight lines to background quasars or galaxies. Halo 
absorption from outflows may be distinguished from gas 
accretion by its location in bipolar flows which emerge 
roughly perpendicular to the galactic disk (Bordoloi et al. 
2011; Bouche et al. 2012; Kacprzak et al. 2012). Second, 
although the emission measure of low density gas typ- 
ically makes halo gas undetectable by direct imaging, 
photons scattered by halo gas can be directly imaged 
(Steidel et al. 2011; Rubin et al. 2011; Prochaska et al. 



2011). Combining the information contained in galaxy 
absorption spectra, resonance emission maps, and sight- 
lines probing galaxy halos should greatly improve the 
accuracy of the estimated mass-loss rates in outflows. 

The most familiar example of resonance scattering in- 
volves the Lya photons produced in H II regions, which 
are scattered in our direction by interstellar and circum- 
galactic gas. P-Cygni Lya line profiles provide direct 
kinematic evidence for scattering by outflowing gas. In 
the simplest model, the redshifted emission component 
is scattered off the back side of an expanding shell while 
the near side imprints blueshifted, resonance absorption 
(Pettini et al. 2002; Verhamme et al. 2008). Since spectra 
of the bluer, less-obscured high-redshift galaxies (Kornei 
et al. 2010) show these P-Cygni line profiles more fre- 
quently than do spectra of more obscured galaxies, the 
detection of diffuse Lya emission in stacks of all types 
of high-redshift, star- forming galaxies was remarkable 
(Steidel et al. 2011). These halos are apparently com- 
mon, extend beyond the UV-continuum emission, and 
indicate that circumgalactic gas scatters Lya photons 
which have escaped from the interstellar medium. 

We suggest a physical analogy exists between the 
prominent Mg II AA2796.35, 2803.53 P-Cygni line profiles 
recently discovered in near-UV spectra of star-forming 
galaxies (Tremonti et al. 2007; Martin & Bouche 2009; 
Weiner et al. 2009; Prochaska et al. 2011; Rubin et al. 
2011; Erb et al. 2012) and the Lya halos detected around 
higher redshift galaxies. The Mg II emission is prominent 
in bluer galaxies and galaxies with relatively low stellar 
mass (Martin et al. 2012), although further analysis in- 
dicates the strongest correlation with specific SFR (Ko- 
rnei et al. 2013). These results appear consistent with 
the ideas that dust attenuation halts the escape of Mg II 
photons and outflows play a significant role in shaping 
the often-observed P-Cygni profiles. 

Motivated by the need to improve empirical mass-loss 
rates and the discovery of Mg II emission around TKRS 
4389 (Rubin et al. 2011, z=0.69), we examined the spa- 
tial extent of Mg II emission in LRIS spectra of indi- 
vidual, blue galaxies at 0.4 < z < lA included in the 
outflow census of Martin et al. (2012) and originally dis- 
covered by the Deep Extragalactic Evolutionary Probe 
2 survey, DEEP2 (Newman et al. 2012). The spectrum 
of the outflow galaxy 32016857 resolves extended Mg II 
emission which we argue most likely reflects stellar con- 
tinuum and H II region emission scattered by the out- 
flow. We also describe extended Mg II emission around 
32010773 and 12019973, two galaxies with spectra indi- 
cating the presence of outflowing gas (blueshifted reso- 
nance absorption and P-Cygni Mg II emission); but, in 
contrast to 32016857, the Mg II surface brightness pro- 
files and kinematics can be fitted with either resonance 
scattering or direct emission from H II regions. Finally, 
the spectrum of 22028686 resolves resonance absorption 
across the galaxy, yet we detect the resonance emission 
over just one side of the galaxy, a situation which may 
provide insight into the conditions required for resonance 
emission to escape. 

We use resonance scattering of the continuum light into 
our sighthne to estimate mass-loss rates for 32016857 and 
TKRS 4389. To outline the physical argument, we con- 
sider spherical outflow geometries subtending a solid an- 
gle Q. Continuum photons emitted at increasingly higher 



frequencies (relative to the resonance) are absorbed by 
atoms at progressively larger radii if the radial veloc- 
ity, v{r), of the outflow increases outwards. At any ra- 
dius, the spatial extent of the scattering region (over a 
fixed frequency range) grows with increases in the ther- 
mal or turbulent velocity. Provided the outflow density 
and the velocity component along the sightline do not 
change much over this interaction region, the computa- 
tion of the scattering optical depth can be simplified us- 
ing the Sobolev approximation (Sobolev & Gaposchkin 
1960; Prochaska et al. 2011). Measurements of the extent 
of the scattering halo and the outflow velocity therefore 
effectively constrain the ionic density of the outflowing 
gas at a radius where the scattering optical depth is close 
to unity. We estimate the ionization correction by calcu- 
lating the photoionization equilibrium and demonstrate 
how the dumpiness of the outflow affects the estimated 
mass-loss rate. 

The paper is organized as follows. We review the data 
in Section 2 and present results for 32016857 in Section 3. 
In Section 4, we examine what galaxy properties favor 
scattered halo gas and present other examples of spa- 
tially extended emission-line gas from our study. We re- 
turn to 32016857 in Section 5, where we use the Sobolev 
approximation to estimate the ionic gas density of the 
outflow, present the photoionization models, and discuss 
the implications for the mass- loss rate. An Appendix il- 
lustrates how the dumpiness of the warm phase of the 
outflow may affect the estimated mass- loss rate. 

We adopt a Hubble constant of Hq = 
70 km s~^ Mpc~^, and we assume a ACDM cos- 
mology with Qm =0.3 and ^^a = 0.7. This cosmology 
yields a scale of 7.884 h:f^ kpc/'' at redshift z = 0.9392. 
We use vacuum wavelengths throughout to refer to both 
ultraviolet and optical transitions and atomic data from 
Morton (2003). Stellar masses and SFR's are derived 
assuming a Chabrier initial mass function (Chabrier 
2003). 

2. DATA 

We previously described 208 Keck LRIS spectra of 
0.4 < z < lA galaxies in Martin et al. (2012). Among 
the 145 spectra with both Mg II coverage and fitted Fe II 
resonance absorption, we have objectively searched for 
Mg II emission (Kornei et al. 2013). We use the parame- 
ters fitted to the Fe II absorption troughs in Martin et al. 
(2012) to describe the shape of the intrinsic Mg II absorp- 
tion trough and definitively detect Mg II AA2796, 2803 
resonance emission in 22 spectra. A higher fraction of the 
spectra may have Mg II emission which is not uniquely 
separated from the absorption troughs. 

The primary spectrum discussed in this paper, 
32016857, was obtained on mask msc32aa (see Table 1 
in (Martin et al. 2012)) using the configuration with the 
D560 dichroic, 600 1 mm~^ blue grism blazed at 4000 
A, and the 600 1 mm~^ red grating blazed at 7500 A. 
Due to the non-photometric observing conditions for this 
mask, absolute flux calibrations for the blue and red spec- 
tra were obtained, respectively, by matching the DEEP2 
photometry in the B and / bands, a procedure which 
also corrects the broad band luminosity for slit losses. 
This absolute flux calibration is less accurate than that 
derived from observations of standard stars under photo- 
metric conditions for most masks; however, all the LRIS 



spectra have accurate relative fluxes independent of ob- 
serving conditions. We refer the interested reader to 
Section 2.2 of Martin et al. (2012) for a more detailed 
description of the instrumental configuration, data ac- 
quisition, and spectral reduction. 

To search for extended emission in the 2D spectral im- 
ages, we fit a Hnear continuum through bandpasses to 
either side of the Mg II doublet. We visually inspected 
the 2D spectra both before and after continuum subtrac- 
tion for any sign of Mg II emission. We collapsed the line 
emission (in both transitions) in the dispersion direction 
to create a surface brightness profile and compared it to 
the instrumental point-spread function (defined by a stel- 
lar spectrum). A similar procedure was followed to create 
an [O II] surface brightness profile from the red spectra. 
Finally, profiles of the continuum surface brightness were 
constructed from both the blue and red spectra directly 
adjacent to the emission line of interest. 

3. CASE STUDY OF THE Z = 0.9392 GALAXY 32016857 

In this section, we focus on the clearest example of 
spatially extended, scattered Mg II emission in our sam- 
ple, 32016857. Figure 1 shows an image of this galaxy 
as well as the blue and red spectra obtained previously 
with LRIS. The i?-band isophotes of 32016857 indicate a 
major axis roughly aligned with the position angle of the 
slit. However, the isophotes clearly show an asymmetry 
along the major axis which is not consistent with the 
projection of an inclined, symmetric disk. Higher resolu- 
tion imaging is required to distinguish whether the ma- 
jor axis reflects the position angle of a disk or some more 
complicated morphology, such as two merging galaxies. 
We wish to emphasize, however, that the surface bright- 
ness distribution and kinematics of the [OH] emission 
are compatible with two emitting regions. In addition, 
we identify the object b" west of 32016857 in Figure la 
as 32016527 in the DEEP2 photometric catalog, which 
suggests the object is a compact galaxy based on its mag- 
nitude and color. 

3.1. Galaxy Properties 

Table 1 lists the measured properties of 32016857. It 
has notably blue U — B color. The SFR estimated from 
the [O II] luminosity, as described in the notes to Table 
1, is approximately 80 M© yr~^. The mass derived from 
the SED fit, log(M/ Mq) = 9.82 (Bundy et al. 2006), 
places 32016857 in the lowest tertile of the LRIS sample 
by stellar mass. This mass estimate should be accurate 
to within a factor of two because the iC-band photome- 
try of 32016857 directly constrains the rest-frame, near- 
IR continuum. At z ^ 1, abundance matching places 
galaxies with stellar mass logM*/ M© = 9.8 in halos of 
total mass logM/,/ M© « 11.5 (Behroozi et al. 2010), 
consistent with the clustering of DEEP2 galaxies which 
places galaxies brighter than Mb < —20.77 in halos more 
massive than logM/,/ M© = 11.3 (Coil et al. 2008). In 
numerical simulations, halos with masses similar to our 
best estimate for 32016857 (specifically Mh < 10^^ Mq), 
gas accretion is dominated by the cold flows never shock- 
heated to the halo virial temperature (e.g., Keres et al. 
2005). 

In the red side spectrum, shown in panel b of Fig- 
ure 1, the resolution of 220 km s~^ FWHM does not 
quite resolve the prominent [O II] AA3727.09, 3729.88 



emission lines. The blue spectrum, 282 km s ^ 
FWHM resolution, reveals strong emission in the Mg II 
AA2796. 35, 2803.53 doublet with a P-Cygni profile (bot- 
tom panels). The recombination lines from hydrogen 
and helium as well as the strong emission in the C II] 
A2324 - 29, C III] A1908.73, [Ne III] A3869.85, A3968.58, 
and [O II] lines come from H II regions photoionized by 
massive stars. The absence of [Ne V] A3426.98 emission 
argues against an AGN as the dominant source of ioniz- 
ing radiation. 

The emission lines in the red spectrum dominate 
the cross-correlation signal with template spectra and 
thereby determine the LRIS redshift of z = 0.9392. The 
root-mean-square (RMS) error in the dispersion solu- 
tion of the red spectrum is only a few km s~^; hence 
the RMS error in the blue dispersion solution, about 
18 km s~^, determines the systematic uncertainty in the 
relative velocities of the low-ionization absorption and 
nebular emission. In the blue spectrum, the centroid of 
the C III] emission line lies just — 30±31 km s~^ from the 
intercombination transition which dominates the emis- 
sion at high density. We readily acknowledge a larger 
discrepancy of 50 km s~^ between the LRIS redshift and 
the DEEP2 survey redshift, likely caused by differences 
in slit position angle but irrelevant to the value of rela- 
tive velocities between the absorption lines and nebular 
emission. 

The absorption lines in the LRIS spectrum are not 
sufficiently resolved to measure the maximum blueshift 
from the Fe II line profiles. The fitted centroid of 
the Fe II absorption lines is blueshifted Vi = — 91 ± 
15 km s~^ (Martin et al. 2012). The projected out- 
flow velocity is more accurately estimated by the two- 
component fit shown in the right column of Figure 2. 
Correcting for the interstellar absorption at the systemic 
velocity increases the blueshift of the Doppler component 
to Vdop = -227 ± 82 km s'^. 

While the A2803 absorption trough may be filled 
in by emission from both Mg II transitions, only 
A2796 emission fills in the intrinsic A2796 absorption 
trough. In Mg II A2803, Figure 1 shows absorption 
blueward of the systemic velocity and redshifted emis- 
sion. The marginally-resolved emission wing extends to 
+310 km s~^ (Martin et al. 2012). The corresponding 
red wing of the Mg II A2796 profile may significantly fill 
in the bluest portion of the intrinsic A2803 absorption 
profile. Emission at the systemic velocity in both transi- 
tions likely fills their intrinsic absorption profiles at the 
systemic velocity. 

Optically thick Mg II absorption will produce an in- 
trinsic A2803 trough with the same equivalent width as 
the A2796 trough which is twice as strong in the opti- 
cally thin limit. Remarkably, however, the absorption 
equivalent width of the Mg II A2796 absorption trough 
is less than that in A2803 as seen in both Figure 1 and 
Table 2. While low S/N ratio could produce such an 
effect through measurement error, our analysis suggests 
a different interpretation. We can understand the net 
Mg II profile provided the intrinsic ratio of emission-line 
strengths Wem(A2796)/Vl^em(A2803) exceed the ratio of 
absorption-line strengths iya6s(A2796)/iya6s(A2803). In 
support of this approach, we note that in Table 2 the net 
emission equivalent width of the stronger Mg II transi- 
tion, —2.13 ± 0.29 A, is higher than in the weaker Mg II 
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Fig. 1. — (a): CFHT i^-band image of 32016857 from Coil et al. (2004). Contours, spaced by a factor of two in surface brightness, 
show an elongated structure; and the slit (PA = 82.0°; 10^0 long by 1^2 wide) runs roughly along the major axis of the galaxy. The 
surface brightness profile along the major axis, however, is clearly more extended to the west than to the east. Due to this asymmetry, an 
inclined disk does not describe the morphology well, and the position angle of the semi-minor axis need not indicate the direction of the 
projected rotation axis, (b): Red spectrum of 32016857 showing strong lines and fitted continuum level. The strong nebular line emission 
determines the systemic velocity, z = 0.939196. Table 2 provides the transition wavelengths, (c): Blue spectrum of 32016857 showing 
strong absorption lines from Fe II, C II] emission, Mg II emission, Fe II* emission, and fitted continuum level. The continuum S/N ratio 
at 2450 A is 8.3 per pixel. Table 2 provides the transition wavelengths, (d): Integrated spectrum of 32016857 zoomed in to illustrate the 
P-Cygni-like line profile of the Mg II doublet. Dotted, vertical lines denote the systemic velocity, and the dashed line marks the normalized 
continuum level. 



line at A2803. Qualitatively at least, we can therefore 
create a line profile with a weak A2796 trough (relative 
to A2803) if the intrinsic absorption troughs are both 
saturated and therefore equal in area. 

To explore this solution quantitatively in a manner 
that fully incorporates the spectral S/N ratio, we fit a 
model consisting of an emission doublet simply added to 
the intrinsic absorption troughs. This model assumes the 
emission region lies beyond most of the absorbing gas. 
The caption of Figure 3 describes the functional form of 
these components. We wrote a custom Markov Chain 
Monte Carlo fitting code in Python to explore the range 
of possible parameter values. The fit shown in Figure 3 
illustrates a typical fit drawn from the Markov chain. 

To estimate the values of the absorption and emis- 
sion components, we used the posterior probability dis- 
tributions compiled for each parameter from the full 
Markov chain to determine the mean values of the 
parameters and their la uncertainties. Our fit de- 
scribes the absorption troughs with a Doppler shift of 
Vabs = — III49 ^^ ^~^' Doppler parameter bjjop = 



81^24 km s ^, covering factor C 



f 



0.9lt^-?^, and op- 



-0.12' 

tical depth at line center tq = 68^44; it finds an emis- 
sion component with a redshift of Vem = lOSl^g ^^ ^~^ 
width a = 14^9^ km s~^, amplitude A2796 
for the A2796 transition, and doublet ratio A2796M2803 
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We find that the intrinsic absorption troughs 
typically have very similar equivalent widths. In con- 
trast, only doublet ratios ^2796/^2803, near 1.5 provide 
acceptable descriptions of the emission. We conclude 
that the optical depth of the emission region in our slit 
is lower than that of the absorbing gas along our sight- 
line. We interpret the difference between the absorption 
and emission doublet ratios as an indication that differ- 
ent physical regions of the galaxy produce the emission 
and absorption - e.g., non-spherical geometries for the 
scattering region and/or significant Mg II emission from 
H II regions. 

The net Fe II absorption troughs more nearly re- 
flect the intrinsic shape of the absorption profile than 
does Mg II absorption. The resonance emission that 
fills in these troughs may be generated by H II re- 
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Fig. 2. — Continuum normalized spectrum showing the near- 
UV, Fe II multiplets UVl, UV2, and UV3. Left column: The 
vertical dashed (dotted) lines mark the Fe II (Fe IP) transitions 
listed in Table 2. The 3a upper limit on the equivalent width 
of the Fe II 2261 hue, not shown, is 0.60 A (Martin et al. 2012). 
Fluorescent emission from the upper level of the A2600 transition 
(z^D^,^) is detected at A2626, and the marginal detections of both 

FeII*A2612 and A2632 mark fluorescent decays from the z^F^^^ 

level, which is excited by absorption in A2587. A net excess of 
emission is also detected near the wavelength of the Fe IPA2396 
transition in UV2 and at FeII*A2365 in UV3. While Table 2 
shows that many of these detections are not individually robust 
statistically, the combination of features clearly argues for strong 
fluorescent emission in the spectrum of 32016857. Right column: 
Two-component joint flt (solid) to the Fe II A 2344, 2374, and 
2587 transitions. The Doppler component (blue, dotted line) has 
a blueshift of Voop = —227 ± 82 km s"-*^ and an equivalent width 
^i:>op(A2374) = 0.41 A. Fit redsiduals (dashed) are shown relative 
to the error spectrum at the bottom of each plot. 



gions or created by the absorption of continuum pho- 
tons in the surrounding gas. As emphasized in Erb 
et al. (2012), H II regions produce much more emis- 
sion in Mg II than in Fe II, so nebular emission will 
have a proportionately larger impact on the shape of 
the net Mg II line profile. In our integrated spectra of 
32016857, the flux ratio of the [O II] and Mg II emis- 
sion, F{[ O II ] AA3726,29)/F(MgII A2796) ^ 42 (un- 
corrected for reddening), is consistent with line ratios 
calculated for nebulae at slightly sub-solar metallicity 
photoionized by a starburst spectral energy distribution. 
For example, comparison of this line ratio to Figure 16 
in Erb et al. (2012), indicates an ionization parameter 
logC/> — 2.9, where the lower limit indicates that the 
reddening-corrected line ratio would require a larger ion- 
ization parameter. In addition, in contrast to the ground 
state of Mg II, which is a singlet, the ground-state of Fe II 
has fine structure. In Fe II, absorption of a resonance 
photon is often followed by the emission of a longer wave- 
length photon (Rubin et al. 2011; Prochaska et al. 2011; 
Martin et al. 2012). This fluorescent emission, clearly 
visible in Figure 2 near the systemic velocity, does not 
fill in the resonance absorption troughs. 

Table 2 summarizes the equivalent widths and veloci- 
ties of Fe II absorption troughs and Fe II* emission. The 
relative strengths of the individual absorption lines are 
not well fit with standard curve-of-growth techniques due 
to differential amounts of emission filling. Using the two- 
component fit shown in Figure 2, Martin et al. (2012) 
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Fig. 3. — Typical example of a fit to the Mg II complex 
in the 32016857 spectrum. We model the emission with two 
Gaussian line profiles sharing a common Doppler shift and ve- 
locity width but independent amplitudes. We add the emission 
component to the fitted absorption doublet described by /(A) = 
1 — Cf -\- Cf exp( — [r2796(A) + T2803(A)]), where the optical depth 
in the bluer line is twice that of the redder line at the same Doppler 
shift; and a Gaussian distribution describes the variation in optical 
depth with Doppler shift, r(A) = tq exp(-[(A- Ao)/(Ao6/c)]2). We 
convolve this joint model with a Gaussian function describing the 
instrumental resolution. The example shown here fit the absorp- 
tion troughs with a Doppler shift of Vabs = —104 km s"-*^, Doppler 
parameter b^top = 70^24 ^m s"-*^, covering factor Cf = 0.97, and 
optical depth at line center tq = 15; it describes the emission lines 
with a Doppler shift of Vabs = 70 km s"-*^, width a = 79 km s"-"^, 
amplitude Ab = 1.67 for the A2796 transition, and doublet ratio 
^b/Ar = 1-67. The doublet ratio fitted to the emission compo- 
nent is larger than that of the absorption component. 

estimate the ionic column density of the Doppler com- 
ponent lies in the range log[NDop{Fe~^)Cf{ cm~^)] = 
14.42 — 15.53, where Cf is a gas covering factor of or- 
der unity. 

3.2. Emission Line Properties along the Slit 

The two dimensional (2D) spectra of 32016857 con- 
strain the spatial extent of the emission regions. Figure 4 
compares the structure of the Mg II and [O II] line emis- 
sion. While the Doppler shift of the [O II] emission shifts 
by approximately 190 km s~^ over 17 kpc along the slit, 
no velocity gradient is detected in Mg II emission. While 
the H II regions producing the [O II] emission proba- 
bly contribute to the total Mg II emission, based on the 
kinematic differences between the 2D spectra, we sug- 
gest that the Mg II emission takes on the velocity of the 
outflowing material, not the velocity of the H II regions, 
because it is scattered by halo gas. 

Figure 5 compares the emission profiles along the slit. 
The continuum emission, rather than showing the sym- 
metric profile of an inclined disk, shows an extended fea- 
ture to the west with strong [O II] but no Mg II emission. 
The [O II] emission is unresolved to the east. In contrast 
to this nebular line emission, the Mg II surface brightness 
profile is spatially extended to the east but unresolved to 
the west. The opposite directions of the extended [O II] 
and Mg II emission along the slit suggest distinct physi- 
cal origins. 

The Mg II surface brightness averaged over a spatial 
resolution element (six pixels) 10.6 kpc eastward of the 
galaxy is 3.1 standard deviations above the background. 
The continuum emission is also spatially extended to the 
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Fig. 4. The 32016857 spectrum near Mg II and [O II]. (a) The Mg II AA2796.35, 2803.53 doublet (vertical hues separated by 770 km s^^) 
is detected in emission and absorption, (b) The [O II] AA 3727.09, 3729.88 doublet (vertical lines separated by 225 km s"-*^) is detected in 
emission, (c) Continuum subtracted Mg II spectrum, (d) Continuum subtracted [O II] spectrum. Linearly spaced contours show a velocity 
gradient along the slit which we attribute to the rotation of the galaxy. 



east. Averaged over the same spatial resolution element 
at 10.6 kpc, the line emission is marginally, 2.1(7, more 
significant than the continuum emission indicating the 
line emission is likely more extended than the continuum 
emission. The IVIg II emission is not spatially resolved to 
the west. 

To estimate the radial extent of the IVIg II emission, 
we introduce a simple model for the surface brightness 
profile measured along the slit inspired by, but different 
in detail from, the approach used by Rubin et al. (2011) 
to model the galaxy TKRS 4389. The emission model 
includes an H II region component and an extended com- 
ponent. The normalized [O II] surface brightness profile 
determines the shape of the H II component. Since the 
observed IMg II profile is spatially resolved to the east 
but not to the west, we require an asymmetric extended 
component. We adopt a semi-circular emission region 
of constant surface brightness. The radius, Rc^ of the 
extended emission describes the projected distance from 
the center of 32016857, and the straight-side of the semi- 
circle is perpendicular to the slit. This 2D model of the 
surface brightness profile is convolved with a Gaussian 
model of the atmospheric seeing, FWHM = 0^8 and 
then "observed" through a V!2 slit. In analogy with Lya 



radiative transfer, photons scattered off the far side of the 
outflow escape as redshifted emission. 

We explore 14 values of Re from 4.9 to 35.9 kpc, cor- 
responding to an angular extent of 0^62 1 to 4V554 along 
the slit. For each model of the extended emission, we fit 
the IVlg II surface brightness profile with a linear combi- 
nation of the H II and extended components. Figure 6 
shows the minimum residuals at R = 11.4 kpc, where 
the separation between models is about 0.34 kpc. The 
extended component and the H II component determine 
the net surface brightness profile to the east and west, 
respectively, in Figure 6. For the best-fit model, the 
eigenvalues of the two components are similar, and the 
extended component contributes 46% of the total IVIg II 
flux. While the residuals show that this model provides 
a reasonable description of the data, it does not tightly 
constrain the radius of the scattering halos. Smaller ha- 
los cannot be excluded with high confidence due to the 
profile smearing caused by atmospheric turbulence, and 
this blurring also means the reduced chi-square values 
increase slowly with increasing halo radius. 

4. GALAXIES WITH HALOS OF SCATTERED EMISSION 

Spectra of about 20% of star-forming galaxies at 2: ~ 
1 show resonance absorption at large (>100 km s~^) 
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Fig. 5. — Surface brightness profiles across 32016857. The line 
emission is normalized to unity. To facilitate comparison to the 
emission line profile, the continuum has been scaled to match the 
maximum line emission. The residuals in the bottom panels show 
the difference of these scaled line and continuum profiles along the 
slit. The point spread function (psf) was measured from a stellar 
spectrum and has been normalized to unity area. Top: The Mg II 
doublet emission is spatially extended to the east, as can be seen 
from the comparison to the stellar profile. At -10.6 kpc, the aver- 
age Mg II surface brightness over one resolution element (6 pixels) 
is significantly (3.1cr) above the background. The continuum emis- 
sion is also extended towards the east. In the lower panel, which 
shows the difference between the normalized line and continuum 
profiles, the line emission is marginally (2.1cr) more extended than 
the continuum emission. Bottom: In contrast, the [O II] surface 
brightness profile is not spatially resolved to the east but is well 
resolved to the west. Since the spatial distribution of the Mg II 
emission does follow that of the [O II] emission, the former clearly 
does not come directly from H II regions. 

blueshifts relative to nebular and stellar lines (IVIartin 
et al. 2012). Since the properties of star- forming galax- 
ies with and without blueshifted absorption do not differ 
very significantly, the absence of blueshifted resonance 
absorption in some spectra has been attributed to the 
collimation of the outflows. A bipolar flow, for example, 
produces a blueshift in the spectrum when the orienta- 
tion of the galactic disk is face-on and the starlight passes 
through the near side of the outflow cone. This interpre- 
tation of the blueshifted fraction suggests outflows are a 
generic property of star-forming galaxies at intermediate 
redshift and raises the question of why scattered emission 
is not detected more frequently. 

Many properties of 32016857 and TKRS 4389 (Rubin 
et al. 2011), the only other distant galaxy detected in ex- 
tended IVIg II emission, are similar. In particular, Rubin 
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Fig. 6. — The normalized Mg II surface brightness profile with 
fitted model (cyan line). This particular model, while not necessar- 
ily unique, provides an acceptable description of the Mg II surface 
brightness with fit statistic x^ = 1.18. The model is a linear combi- 
nation of the emission from an HII region component (red, dotted 
line) and an extended halo (blue, dashed line). The normalized 
[O II] surface brightness profile defines the HII region component. 
The extended component is modeled by a semi-circular region of 
radius Rq and constant surface brightness. The best-fit radius for 
this scattered component is Rq = 11.4 kpc, and the spacing be- 
tween models is 0.3 kpc. The lower panel shows that the errors in 
the data values are comparable in magnitude to the the residuals 
from the fit (data minus model). The inset, however, shows the re- 
duced chi-square values over a large range in Rq- Since the range 
in Re over which x^ increases by 1 from its minimum value at 
11.4 kpc is quite large, models with halo radii considerably smaller 
or larger than 11.4 kpc are not ruled out. 

et al. (2010) suggest that TKRS 4389 is exceptionally 
bright for its redshift oi z = 0.6943. It is slightly bluer 
and more luminous than 32016857 and has about three 
times as much stellar mass. A rough estimate for the 
SFR is ^ 80 M© yr-^ The deepest parts of the Fe II ab- 
sorption troughs are blueshifted roughly 200 km s~^. An 
HST image of TKRS 4389 shows signs of recent merger 
activity. Rubin et al. (2010) cite the weak [Ne V] emis- 
sion, Wr = -0.84 ± 0.1 A, as evidence that TKRS 4389, 
in contrast to 32016857, hosts a low-luminosity active 
galactic nucleus (AGN). 

The very blue colors of 32016857 and TKRS 4389 are 
typical of Mg II emitters, which have lower masses, bluer 
colors (i.e., less reddening), and higher sSFR than the av- 
erage star-forming galaxy at intermediate redshifts (Erb 
et al. 2012; Martin et al. 2012, Rubin et al., in prep). 
However, these two galaxies are more luminous than the 
average Mg II emitter; their blue luminosities fall in the 
highest tertile of the Martin et al. (2012) sample. In 
spite of the high luminosity of 32016857, its stellar mass 
falls in the lowest tertile of our sample. The mass of 
TKRS 4389 is significantly higher, near the upper bound 
of our middle tertile. Galaxies that are both as lumi- 
nous and as blue as 32016857 are rare. Only 9 of 208 
galaxies in the parent sample have U — B < 0.459 and 
Mb — 5 log /i 70 < —21.3. Three of these 9 galaxies are 
Mg II emitters by the definition of Kornei et al. (2013), 
but only the 32016857 spectrum shows spatially extended 
Mg II emission. 

Many factors, including the size of the spectroscopic 
aperture, the solid angle of the outfiow and spectroscopic 
aperture, the amount of interstellar gas and dust, and the 
orientation of the galaxy and wind affect whether an out- 



flow can be detected in resonance emission. Accordingly, 
it is unclear whether the resonance photons only escape 
from these galaxies or whether the emission is otherwise 
spread over too large a solid angle to detect it through 
a narrow slit. Drawing from our observation that both 
32016857 and TKRS 4389 have uncommonly large lumi- 
nosities, mapping Mg II emission in and around galaxies 
like these nine may be the most direct way to determine 
whether scattered halo emission is a generic property of 
intermediate redshift galaxies. While this could be ac- 
compHshed with LRIS using additional longslit position 
angles, the next generation of optical spectrographs with 
large integral field units, particularly KCWI at Keck and 
MUSE at VLT, will be well positioned to characterize the 
halo emission. 

With the aim of better informing these future investi- 
gations, we describe in Section 4.1 additional highlights 
from our investigation of the spatial extent of Mg II emis- 
sion in deep, LRIS spectra. From the union of a color 
selected sub-sample and the sample selected by Mg II 
emission in the integrated spectra, we identified two more 
(for a total of three) examples of extended Mg II emis- 
sion which we discuss in this section. We also describe a 
noteworthy example of extended Mg II absorption with 
unresolved emission in 22028686. Finally, in Section 4.2, 
we remark on the absence of Fe II A2383 emission. We 
return to 32016857 in Section 5 where we discuss the im- 
plications of scattered Mg II emission for mass-loss rates. 
The reader primarily interested in the physical proper- 
ties of outfiows rather than the challenge of detecting 
scattered emission should proceed directly to Section 5. 

4.1. Additional Examples of Spatially Resolved Mg II 
Structure 

The bluest fifth of the Martin et al. (2012) sample con- 
tains 42 galaxies with U — B < 0.459 and is of particular 
interest since the Mg II emitters tend to have blue col- 
ors (Martin et al. 2012; Kornei et al. 2013). Ah but 4 of 
these spectra cover the Mg II doublet, yet only 10 of these 
38 show prominent Mg II emission in integrated spectra. 
Among this subsample, we found only two examples of 
spatially extended emission in the Mg II doublet. In ad- 
dition to 32016857, the Mg II emission in our spectrum 
of 32010773 is spatially extended. 

Figure 7 shows the surface brightness profiles of the 
Mg II emission (panel d), blue continuum, [O II] emis- 
sion (panel c), and red continuum along the slit cross- 
ing 32010773. In contrast to the surface brightness pro- 
files measured across 32016857, the Mg II emission across 
32010773 is no more extended than the blue continuum 
emission. In the continuum-subtracted, 2D spectra, the 
velocity gradients of the Mg II and [O II] emission along 
the slit do not obviously differ, so the Mg II emission may 
come directly from H II regions. Our data do not require 
(nor do they rule out) scattered Mg II emission. This 
galaxy does clearly have an outfiow; the large blueshift 



of the Fe II resonance lines, Vi 



-198 ±46 kms-\ is 



quite significant. The color of 32010773, U - B = 0.23 
is among the bluest in the sample; but the luminosity of 
32010773, Mb - 51og/i70 = -20.68, is much lower than 
the luminosity of either 32016857 or TKRS 4389. The 
absence of resonance emission that is clearly extended 
relative to the continuum around 32010773 is consistent 
with the conjecture that detection of a scattering halo 



requires especially luminous galaxies. 

In addition to these 10 Mg II emitters with very blue 
color, Kornei et al. (2013) identified 12 more Mg II emit- 
ters in the Martin et al. (2012) sample with U — B red- 
der than 0.459. Among these redder Mg II emitters, we 
found spatially extended Mg II and continuum emission 
across 12019973. Figure 8 shows the variation in the 
Mg II and [O II] emission along the slit. Across 12019973, 
much like 32010773, the surface brightness profile of the 
line emission is no more extended than the continuum 
emission, neither proving nor ruling out scattering by an 
outfiow. We note that the Doppler shift of the Fe II res- 
onance lines, Vi = —59 ± 28 km s~^, provides marginal 
evidence for an outfiow from this galaxy. The color, 
U - B = 0.539, and stellar mass, log(M/ M©) = 10.15, 
of 12019973 are typical of the sample. This galaxy has a 
relatively low luminosity, M^ — 5 log hro = —20.57, which 
may contribute to the faint ness of any extended halo. 

Figure 9 shows the 2D spectra of another Mg II emit- 
ter, 22028686, a star-forming galaxy with U-B > 0.459. 
Although the Mg II emission is not spatially resolved, 
the resonance absorption troughs extend well beyond the 
emission region. The Mg II absorption lines are tilted in 
the opposite sense of the velocity gradient indicated by 
the [O II] Doppler shifts along the slit. Since the [O II] 
emission defines the systemic redshift, the redshifted and 
blueshifted [O II] emission to the northeast and south- 
west of the galactic center reveal either the projected 
rotation of the galaxy or the orbital motion of two merg- 
ing galaxies. The centroid of the Mg II absorption to 
the southwest of the galaxy is consistent with the veloc- 
ity of the blueshifted side of the galaxy. It is the even 
larger blueshift of the Mg II absorption centroid to the 
northeast that makes the absorption lines tilt in the op- 
posite direction of the [O II] emission. Where the Mg II 
absorption is blueshifted on the northeastern side of the 
galaxy, we detect strong resonance emission at redshifted 
velocities. No resonance emission is detected along the 
southwestern half of the slit. We note that the Fe II ab- 
sorption in the integrated spectrum does not show a net 
Doppler shift, Vi = -1.9 ±9.2 km s'^. Much of the shift 
in the centroid of the Mg II absorption may therefore be 
caused by variations in emission filling along the slit, and 
the spatial gradient in the absorption velocity likely does 
not refiect galactic rotation. 

The detection of Mg II emission from just a portion 
of the absorbing region in 22028686 may help illuminate 
the physical conditions required for the resonance emis- 
sion to escape the galaxy. To gain this insight, we re- 
view studies which have mapped Na I resonance absorp- 
tion across galaxies. As illustrated for example by spec- 
tral mapping of NGC 1808 (Phillips 1993) local ULIRGs 
(Martin 2006), and Mrk 231 (Rupke et al. 2005; Rupke 
& Veilleux 2011), the absorption from outfiows is not 
redshifted anywhere across the galaxy because the ab- 
sorbing gas must lie in front of the continuum source. 
Across local ULIRGs, the Na I velocity gradient is very 
similar to the slope of the rotation curve, but the ab- 
sorption trough is highly blueshifted everywhere along 
the slit. This situation would be expected if the outfiow 
was launched from a rotating disk, and we might expect 
a small velocity gradient in the absorption Doppler shift 
(not seen) comparable to the tilt of the nebular emission 
lines. In spite of the very limited spatial resolution across 
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Fig. 7. Observations of 32010773. (a) The galaxy, shown in a CFHT R-band image, was observed with LRIS through a V/2 by 10(^0 
slitlet at a position angle of 82.0° (on mask msc32aa as described in Table 1 of Martin et al. (2012)). (b) The blue LRIS spectrum reveals 
a P-Cygni Mg II profile, (c) Normalized surface brightness profiles of the [O II] emission the nearby continuum are spatially extended 
relative to the point-spread function, (d) Normalized surface brightness profiles of the Mg II emission is partially resolved but no more 
spatially extended than the blue continuum emission. The Gaussian FWHM fitted to the Mg II surface brightness profile, 1(^25 ±0^11, is 
not significantly larger than the FWHM of the blue continuum profile, 1(^14 zb 0(^02. The fitted centroid of the Mg II emission profile is 1.2 
pixels (0(^16) east of the continuum centroid, a 2.9a discrepancy in position. 



22028686, we speculate that the geometry might instead 
mimic NGC 1808 where differences in extinction shape 
the resonance hne profile of the outflow. As illustrated in 
Figures 1 and 4 of Phillips (1993), the blueshifted lobe of 
the outflow is detected in absorption across the inclined 
disk. Resonance emission is detected from the far side 
of the outflow where it pokes out from behind the disk 
but plumes of dust visible on the near side of the outflow 
at the location of the slit presumably destroy the reso- 
nance photons. Hence the implication from the 22028686 
spectrum seems to be that the sightline towards the ab- 
sorbing gas on the northeast side encounters less dust 
than that towards the southwest side of the slit. Higher 
resolution imaging could confirm the presence of a disk, 
determine the disk orientation, and perhaps reveal the 
locations of dust filaments in this galaxy. 

In summary, among 145 spectra with coverage of the 
Mg II doublet, including 50 galaxies with either promi- 
nent Mg II emission or very blue U — B color, we find 
only three examples of spatially resolved Mg II emis- 
sion: 32016857, 32010773, and 12019973. Our analysis 
of 32016857 doubles the number of known galaxies with 



Mg II emission extending beyond the stellar continuum. 
With only two examples where the extended emission 
must be scattered, we can only begin to speculate about 
the conditions producing scattered Mg II halo emission. 
However, since the unique properties of 32016857 and 
TKRS 4389 are their high luminosities (for their blue 
color), we suggest that large luminosities are essential 
to raise the surface brightness of the scattered emission 
to detectable levels. Further support for the idea that 
most scattering halos are below current detection lim- 
its comes from stacking the surface brightness profiles of 
1 < z < 2 Mg II- emitting galaxies; their Mg II emission 
is marginally more extended than the near-UV contin- 
uum at radii of 0^8 (Erb et al. 2012). 

4.2. Absence of Extended Emission in Fe II Lines 

Our spectrum of 32016857 shows prominent fluorescent 
emission in Figures 1 and 2. We did not find extended 
emission in resonance lines of Fe II or fluorescent Fe II* 
lines nor did Rubin et al. (2011) in their TKRS 4389 
spectrum. The Fe II A2383 absorption is optically thick 
in most galaxy spectra; under LS coupling, the only 
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Fig. 8. — Observations of 12019973. (a) The galaxy is shown on the V band HST image described in Kornei et al. (2012). We obtained 
the LRIS spectrum through a 1(^2 by 10(^0 slitlet at a position angle of 170.0° (on mask mscl2ee as described in Table 1 of Martin et al. 
(2012).) (b) The blue LRIS spectrum reveals a P-Cygni Mg II profile, (c) Along the slit, the strong [O II] emission is slightly more extended 
than our model of the psf but not significantly more extended than the red continuum emission, (d) The Mg II surface brightness profile 
exemplifies the low S/N ratio typical of Mg II emission. 



permitted decay from the upper energy level, here the 
excited z^F^-^,^ state of Fe+, is to the ground a^Dg^2 
state. The absence of Fe II A2383 emission in spectra 
with prominent Mg II A2803 emission is at first surprising 
considering the similar oscillator strengths of the A2383 
and A2796 transitions and the (only) 20% larger cosmic 
abundance of magnesium (relative to iron) . 

Several factors may contribute to the paucity of Fe II 
A2383 emission. A greater depletion of iron (relative to 
magnesium) onto grains, for example, may reduce the 
Fe II A2383 optical depth relative Mg II A2803. In the 
Milky Way, however, the relative depletion of Fe and Mg 
in the halo clouds is much lower than the ratio mea- 
sured for disk clouds (Savage & Sembach 1996). Hence, 
we might anticipate a similar depletion of Fe and Mg 
onto grains in outflows. Alternatively, since iron is not 
a light element, deviations from pure LS coupling might 
in principle allow decay to another term thereby reduc- 
ing the resonance emission (Shull, pvt. comm.), but we 
find no observational evidence for unidentified emission 
lines in support of this conjecture in the near-UV spec- 
tra. Third, as suggested recently by Erb et al. (2012), 
boosting the Mg II emission, rather than suppressing the 
Fe II emission, offers another explanation. Those authors 



show that photoionized gas produces more emission in 
Mg II than Fe II, so emission from H II regions may boost 
the intensity of the scattered Mg II radiation relative to 
the scattered Fe II photons. 

The most important factor driving this discrepancy, 
however, appears to be the ionization corrections. Gi- 
avalisco et al. (2011) show (see their Figure 14) that pho- 
toionization models of very low column density {Nh ~ 
10^^ cm~^) gas indicate a large plausible range in the 
relative column density; Nj^gii/Npeii increases from 
approximately unity {logU = —4) to just over 100 at 
\ogU = -2. Even though the first (7.646 eV, 7.870 eV) 
and second (15.035 eV, 16.18 eV) ionization potentials 
of Mg and Fe, respectively, are quite similar, the Fe is 
mainly in Fe^^ while the Mg is in Mg^^ in our pho- 
toionization models, which we describe in the next sec- 
tion. This happens because it takes only 30.6 eV to ionize 
Fe++ to Fe+^ (and 54.8 eV to make Fe^"^) but 80 eV to 
remove another electron from Mg^^. 

5. DISCUSSION 

In this section we apply the idea of resonance scatter- 
ing to the spatial extent of Mg II emission. We take as a 
model for the outflow a purely radial flow as illustrated in 
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Fig. 9. Observations of 22028686. (a) CFHT R-band image. The galaxy was observed with LRIS through a V/2 by 15^0 shtlet at 
a position angle of 25.0° on mask msc22bb as described in Table 1 of Martin et al. (2012). (b) A 140 A wide segment of the blue LRIS 
spectrum centered on the Mg II doublet. Wavelength increases to the right, and northeast is up. The prominent absorption lines exhibit 
a velocity gradient across the galaxy; the centroid shifts blueward from the systemic velocity on the southwest side to approximately 
— 120 ±30 km s"-*^ 1(^08 (8.31 kpc) to the northeast, (c) Continuum subtracted blue spectrum. The Mg II emission component is spatially 
unresolved. The continuum emission subtends a larger angle than the resonance emission, (d) Integrated Mg II profile. The P-Cygni 
shape of the line profile, blueshifted absorption and redshifted emission, is typical of our Mg II emitter sample, (e) A 90 A wide segment 
of the red LRIS spectrum. The blended [O II] doublet is spatially extended. Wavelength increases to the right and northeast is up. 
Linearly spaced contours show a slight tilt of the nebular emission in the opposite sense of the velocity gradient apparent in panel (b). (f) 
Continuum subtracted red spectrum. The faint, serendipitous emission line near the northeastern end of the slit is 4^44 northeast of the 
peak continuum emission along the slit, and we identify it as [O II] emission from the R = 24.8 DEEP2 galaxy 22028464, which is visible 
4(^66 northeast of 22028686 in panel (a). 



Figure 10. The radial component of the velocity is crit- 
ical since the photons emitted by stars and H II regions 
emanate from smaller radii. Atoms in the circumgalac- 
tic medium (CGIVI) scatter photons at the resonance fre- 
quency in their rest frame. Outflowing gas absorbs pho- 
tons blueward of resonance. The Doppler shift of the 
absorption trough directly reflects the velocity of the gas 
located between the galaxy and the observer. The red- 
shifted emission from the far side of the outflow may es- 
cape and produce a redshifted emission line. The size of 
the scattering halo indicates the radius where the optical 
depth along a radial ray is approximately unity. 

The clouds accelerate near the galaxy due to the ra- 
dial decline in both thermal gas pressure and radiation 
pressure. The continuity equation gives the mass flux, 

M{r)=nfcp{r)v{ry, (1) 

at radius r through the solid angle ft subtended by the 
wind. The cloud covering fraction, /c, describes the frac- 
tion of the outflow solid angle covered by low-ionization 
clouds (as seen from the galaxy). 

Defining ffiion as the mass in the warm outflow per ion, 
the column density between the observer and the launch 



radius, Rq^ is 



N,, 



p{r)/miondr. 



(2) 



Ro 



Solving Eqn. 1 for /)(r), we see that the outflow density 
decreases with the inverse-square of the radius for the 
simple case of a constant velocity outflow with no mass- 
loading beyond the launch radius. For this example, the 
integral in Eqn. 2 yields a column density that declines 
inversely with the launch radius Rq. Accurately mea- 
suring the ionic column density will still leave the mass- 
loss rate uncertain at the level of the plausible range of 
launch radii as well as the total mass per ion. Increases in 
column density strengthen the absorption trough signifi- 
cantly in the optically thin limit; but for saturated lines 
like the Mg II doublet, the range in gas velocity along the 
sightline largely determines the equivalent width of the 
intrinsic absorption troughs. The measured absorption 
equivalent width therefore provides only a lower limit 
on the ionic column density. The absorption troughs 
imprinted by outflows in galaxy spectra provide no con- 
straint on the location of the gas along the sightline. 
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Fig. 10. — Schematic diagram of a spherical outflow around a galaxy illustrates an opening angle of ^2 = 47r steradians, equivalent to a 
conical outflow with Ocone = tt. The regions colored red and blue denote redshifted and blueshifted emission, respectively; and the hashed 
region denotes the location of atoms producing blueshifted absorption. For a delta- function line proflle at the resonance uq^ atoms at each 
point along a sight line through the outflow absorb continuum photons at a unique frequency u = i/o(l + v{r) / c)~^ ^ where r = \/\P' + z^ 
along a sight line at impact parameter h. The emitted line photons have Doppler shifts corresponding to the range of projected outflow 
velocities, ^'(r)cos^. For a constant velocity gradient, the density falls as r~^, and we flnd the maximum Sobolev optical depth along the 
sightline at z = 0. 

Resolving scattered line emission constrains the loca- 
tion and density of the gas flow. The scattering atoms 
must lie at radii at least as large as the projected impact 
parameter. The scattering optical depth scales with the 
gas density for a given model (or measurement) of the 
velocity gradient in the outflow. 

In this section, we first we estimate the ionic density 
and radius at the scattering surface in 32016857 and 
TKRS 4389. Then we discuss the limits on the relation- 
ship between the total gas density and the ionic density. 
Finally, we illustrate the implications for the mass flux 
in low- ionization gas. 

5.1. Density of Mg^ in Outflowing Gas 

We model the outflows with a radial velocity field. For 
32016857, we picture a scenario where our slit intersects 
the outflow cone on only the west side of galaxy. In 
TKRS 4389, in contrast, the longslit observation resolves 
Mg II emission on both sides of the slit (Rubin et al. 
2011). A simple model with a constant surface bright- 



ness halo suggests the impact parameter of the scatter- 
ing surface is 8.25 < 6(kpc) < 13 for TKRS 4389 (Rubin 
et al. 2011) and roughly 11.4 kpc from 32016857. 

The interaction region where a resonance transition 
such as Mg II A2796 can be scattered depends on the 
line profile. To simplify the calculation, we consider a 
delta function line profile centered on the resonance fre- 
quency uq (in the reference frame of the bulk flow). Pho- 
tons emitted by the starburst at frequency v propagate 
radially until they encounter a gas parcel with radial 
velocity v{rs) = c{iy — i^o)/^o- We call the interaction 
spot the Sobolev point, rs- Provided velocity v^ /i, and 
the radial velocity gradient change slowly around r^, the 
Sobolev expression (Sobolev & Gaposchkin 1960; Lamers 
& Cassinelli 1999) describes the scattering optical depth. 
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(3) 



where / represents the transition oscillator strength.^ 

Note that the stimulated emission term has been dropped from 
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The incident photons and the velocity of the bulk flow 
are both radial, so the cosine of the angle 6 between them 
takes the value ji = cos^ = 1. The frequency of the pho- 
tons absorbed from the starburst varies only with v{r) 
for a spherical outflow. 

With /i = 1 in Eqn. 3, the Sobolev optical depth for 
the stronger Mg II transition is given by 



'2796 



(r) = 4.57 X 10" 



-7 '^ - 

cm s 



^Mg- 



.(r) 



(4) 



At some radius along the sightline at impact parameter 
h ^ 11 A kpc, we expect the scattering optical depth for 
the starburst photons to be approximately unity. We es- 
timate the gas density by setting T^96(r) ~ 1 and obtain 



^Mg+ ^ 1-4 X 10 



-9 cm-3 



dv 



11 A kpc 



230 km s" 



dr 



.(5) 



If the sensitivity of our spectrum, rather than the ra- 
dius of the last scattering surface, determines the ob- 
served extent of the scattered emission along the slit, 
then the optical depth (and the inferred density) would 
be higher than this estimate. However, we think mod- 
ehng the spatial extent with the last scattering surface 
is roughly correct because the emission is not optically 
thick at 6 ~ 11.4 kpc; at the largest impact parameters 
where Mg II emission is detected, only A2796 is detected, 
consistent with a two-to-one flux ratio and an optical 
depth in A2803 less than unity. The gas density esti- 
mated from Eqn. 5 depends on the velocity gradient at 
the Sobolev point, and this radius is not uniquely deter- 
mined from the impact parameter h alone. 

Fig. 10 illustrates the geometry for the outgoing, scat- 
tered photons at impact parameter h. Along a sightline, 
the angle ranges from 90° where the flow is perpendic- 
ular to the sightline to a minimum, Omim and maximum 
Omax = 180° — Omim Constrained by the Doppler shifts 
of the absorption and emission lines. The Doppler shift 
of a scattered photon reflects the projected velocity com- 
ponent, Vz^ and will be most blueshifted (redshifted) as 
Omin -^ 0° {Omax ^ 180°). Thcsc cxtrcma mark the 
scattering that occurs at the largest radius along a sight- 
line; and we infer the largest (smallest) radius from the 
lower (upper) limits on Omin- The error bars on the angle 
come from considering the projection of the radial ve- 
locity (measured in absorption) onto to the line-of-sight 
component measured at impact parameter b in emission. 
For 32016857, we estimate 38° < Omin < 75°, where 
the spectral fitting constrains the most blueshifted ab- 
sorption ((|VdopI + b)) to the range 185 — 383 km s""*^; 
and the the redshifted emission (Vz + a) extends to 
at least 100 — 145 km s~^. The maximum radius from 
which we detect scattered emission is therefore between 
1.046 ^ 12 kpc and 1.616 ^ 18 kpc. 

We estimate an average velocity gradient of 
230 km s""*^ over 12 to 18 kpc. Substituting these 
values into Eqn. 5, we infer an ionic density between 
8.9 X 10-^^ cm-3 for rs = 18 kpc {Omin = 38°) up 
to 1.3 X 10-^ cm-3 for rs = 12 kpc {Omin = 75°). 
A specific form for the acceleration of the outflow, 
v{r) = vo^/\nr/Ro, demonstrates how a more phys- 
ical model might change our density estimate. The 

Eqn. 3 due to the low gas densities in galactic outflows. 



acceleration of a momentum driven outflow begins 
rapidly and then slows down with increasing distance. 
This model reaches a characteristic velocity vq at a 
radius a few times the launch radius, r = 2.72Rq. 
Choosing Rq such that v{rs) = '^o and rs = 2.72i?o, 
we obtain dv/dr = O.bvo /rs- Placing the Sobolev point 
at 2.72Rq for comparative purposes, the estimated 
velocity gradient at rs takes a value half as large as the 
constant acceleration limit; and our density estimate 
would decrease by a factor of two. 

For TKRS 4389, Rubin et al. (2011) measured the 
Doppler shift of the low-ionization absorption lines and 
fit a radius for the scattering halo to the surface bright- 
ness profile of the Mg II emission. Substituting their val- 
ues of Vdop ^ 200 - 300 km s"^ and b ^ 8.25 - 12.4 kpc 
into Eqn. 5, we estimate nMg-\- ^ 1.1 — 2.6 x 10~^ cm~^. 
Since the Sobolev radius could lie at radii as large as 
rs ^ 1.56, where we have used the Doppler shift of 
the A2803 emission, Vz ^ 226 km s~^ in Table 1 of 
Rubin et al. (2011), and v = 300 km s~^ to estimate 
Omin = cos~^{vz{rs)/v{rs)) ^ 41°, we estimate a lower 
limit on the ionic density of roughly 0.73 x 10~^ cm~^. 
The slightly larger density, nMg-\- ^ 7x 10~^ cm~^, found 
by Rubin et al. (2011) is consistent with substituting the 
-800 km s~^ Doppler shift of the blue wing on the Mg II 
profile into Eqn. 5 along with the minimum Sobolev ra- 
dius. These arguments constrain the ionic density at the 
Sobolev radius in the TKRS 4389 outflow to the range 
0.7- 7 X 10-^ cm-3. 

5.2. Gas Density 

Given the ionic density nMg+, the total gas density 
scales inversely with the ionization fraction, metallicity, 
and dust depletion. Defining the ionization fraction 



x(Mg+) =nMg+/nMg, 
the hydrogen number density is 



riH 



7?(Mg)d(Mg)x(Mg+)' 



(6) 



(7) 



We scale to solar metallicity, 7^(Mg) = 3.8 x 10~^, which 
is appropriate (to a factor of two or better) for a wind 
loaded with interstellar gas from a galaxy with the stel- 
lar mass of the 32016857 (Zahid et al. 2011). We chose a 
fiducial value for the depletion of Mg onto grains typical 
of clouds in the Milky Way disk, c^(Mg) ^ 6.3 x 10"^ 
(Savage & Sembach 1996), and note that a lower deple- 
tion similar to clouds in the Milky Way halo {d ~ 0.28) 
would lower our density estimate by a factor of four. 

To estimate x(Mg+), we start with a rough assump- 
tion, which we relax later, namely that the bulk of the 
Mg ions are either Mg^ or Mg^^. Then 



x(Mg+) ^riMg+friM^ 



.. = ^^^&^U-\ (8) 



CCr(Mg+) 

where the ionization parameter U is defined as 

Q 



u 



Airr'^crip ' 



and 



^Mg+ ~ ^ 4. /^ X iU . 

CCFMg+ 



(9) 
(10) 
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In the last expression, the recombination rate 
aMg+{T) ^ 3.5 X 10-12 cm^ s'^ at T = 10"^ K 
(Shuh & van Steenberg 1982), and the photoionization 
cross section (JMg-\- ^ 2.5 x 10"^^ cm^ (Verner et al. 
1996). For 32016857, we estimate the ionizing photon 
luminosity Q ~ 1.3 x 10^^ s~^ from the Ha luminosity 
(see description in Table 1) and the relationship in 
Kennicutt (1998). 
Using Eqn. 8, we can solve for uh ^ ^e 



riH 



f^Ug+Q 



1/2 



(11) 



At the location where r2796 
a H density 



47rceMg+/?(Mg+),d(Mg)r2_ 

1, the fiducial values give 



uh = 0.30 cm 

6.3x10-2^'/' 



( ^Mg+ \ 

vio-9 cm-3y 



1/2 



10 
3.8 X 10-^ 



11.4 kpc 



rs 



1/2 



Q 



1055 g-1 



1/2 



(12) 



d{Mg) J V ^(Mg) 

The implied ionization parameter at the Sobolev point 
is 

'0.30 cm-^^ 



U{rs) ^ 7.6 X 10" 



1055 g-1 



Tie 



11 A kpc 



rs 



.(13) 



The resulting ionization fraction of Mg^, 



x(Mg+) ^ 1.4 X 10" 



10' 



55 



( ^- ) 

V0.3 cm-3y 

rs 



11.4 kpc 



(14) 



indicates most of the Mg in the warm phase of this out- 
flow has been ionized to Mg^^. For comparison, inspec- 
tion of Fig. 5 in Murray et al. (2007) shows that x(Mg^) 
lies between 0.1 and unity over broad range in gas den- 
sity in ULIRG outflows due to the very low ionization 
parameter in dusty galaxies. 

Using Cloudy version 13.00 (Ferland et al. 2013), we 
calculated the photoionization equilibrium of a constant 
density gas slab illuminated by a starburst and examined 
x(Mg^) as a function of U. Over the range 8 x 10~^ < 
U < 0.08, a simple power law flt gives 



x(Mg+) ^ riMg+friMg = eMg+U~ 



(15) 



with eMg-\- = 4 X 10 ^ and a = 0.94. Dropping the 
assumption introduced in Eqn. 8, Equation 11 is modified 
to 



uh 



^Me+ 



7?(Mg+)d(Mg)eM^+ 



l/(l+a) r 



47r cr^ 



a/(l+a) 



(16) 



For these improved values of x(Mg^) and a, Eqn. 12 
becomes 



riH = 0.15 cm ^ 

6.3 X 10-^ 
^(Mg) 



( ^Mg+ \ 

vio-9 cm-3y 



0.52 



11.4 kpc 



rs 



0.96 



0.52 



3.8 X 10" 

^(Mg) 



0.52 



Q 



10^5 s- 



0.48 

(17) 



In the previous section, we argued that a sightline at 
impact parameter b ^ 11.4 kpc through the outfiow in- 
cludes emission from rs ^ b up to rs ^12 — 18 kpc. For 
this range of radii and nMg-\- = 1.4 x 10~^ cm~^, our best 
estimate for the hydrogen density at the Sobolev radius 
becomes 



UH ^0.12-0.18 cm" 



(18) 



To facilitate comparison to TKRS 4389, we have 
adopted the same fiducial corrections for abundance and 
depletion as Rubin et al. (2011). In contrast, however, 
application of our argument for the ionization correc- 
tion to TKRS 4389 will yield a significantly larger H 
density than the Rubin et al. (2011) estimate of uh ^ 
0.003 cm~^ because they assumed x(Mg^) ^ 1. Us- 
ing the ionic density nMg+ = 2.6 x 10~^ cm~^ (from 
the previous section) in Eqn. 12, we estimate uh = 
0.28 - 0.67 cm-3 for r^ = 20 - 8 kpc. Over this same 
range of plausible values for r^, our more accurate es- 
timate from Eqn. 17 is uh = 0.15 - 0.33 cm~^. We 
conclude that the gas density at the Sobolev points in 
the TKRS 4389 and 32016857 outfiows are very similar. 

5.3. Mass-loss Rate 

Using the H gas density from the previous section and 
assuming 1.4 atomic mass units per hydrogen atom, the 
mass loss rate in the warm gas from Eqn. 1 becomes 

*W = ^""-^'-(")(f)(TTfFpc,' 

l230kms-V U.3cm-») ' *''' 

We have used the observation that resonance lines are 
blueshifted 100 km s~^or more in 20% of star- forming 
galaxies at z ~ 1 (Martin et al. 2012), and the in- 
terpretation that the outfiows must subtend roughly 
0.2 x47r = 0.87r ^ tt steradians on average. For 32016857, 
the mass fiux in the (smooth) warm outfiow is therefore 
approximately 330 - 500 M© yr~^, which is 4 - 6 times 
larger than the SFR of 80 M© yr"^ (for the Chabrier 
IMF). For TKRS 4389, we estimate a mass fiux of 35 to 
40 M© yr-^ 

We can use the empirical limits on the Mg II column 
density as a consistency check on the above argument. 
As can be seen in Figure 3, resonant emission fills in the 
Mg II absorption troughs, so the absorption equivalent 
width must be estimated from a fitted model. Even the 
weaker transition of the Mg II doublet is optically thick 
in any reasonable fit to the 32016857 spectrum; and, in 
Figure 3, it has an equivalent width of VF2803 = 2.04 A. 
We can use the linear relationship between W2803 and 
column density, valid in the optically thin limit, to place 
a lower limit on the ionic column density of A/'(Mg^) > 
1.0 X 10^"^ cm-2. 

Substituting Eqn. 1 into Equation 2 defines the rela- 
tionship between this column density and the mass-loss 
rate. For a constant velocity outfiow, v{r) = Voop^ it 
yields 

M ^ ftfhionfcSNionVDopRO' (20) 

The mass-loss rate estimated from absorption lines there- 
fore scales linearly with the inner radius of the outfiow. 
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Rq, which physicahy represents a launch radius. Using 
our mass-loss rate estimate from the scattered emission 
and the lower limit on the column density, however, we 
obtain an upper limit on the maximum launch radius 



Ro < 2.0 kpc 



M 



7T\ fix 10-^^ g 



500 Mq yr" 



10^ 




m 



230 km s- 



NMg+ J V yOop J 

Dynamical models of outflows relate the launch radius to 
the disk scaleheight (Mac Low et al. 1989; Murray et al. 
2011), suggesting values of Rq ranging from ^ 100 pc 
for thin gas disks to several kiloparsecs for dwarf galax- 
ies and very turbulent systems. We conclude that the 
strength of the saturated absorption lines, while in no 
way providing a test of our large ionization correction, 
are fully consistent with the ionic density of Mg^ in- 
ferred from the scattered emission. 

Up to this point, we have tacitly assumed a smooth 
outflow with /c ^ 1. Higher resolution spectroscopy 
of outflows, however, indicates that saturated, low- 
ionization lines are not always black; and this partial 
covering of the continuum light indicates the warm out- 
flows do not always have unity covering fraction (Mar- 
tin & Bouche 2009). From a physical standpoint, inter- 
stellar gas is likely entrained in a hotter wind to create 
the outflow. In the Appendix to this paper, we outline 
simple physical arguments for the density and number 
of warm clouds in outflows. We find that the clumps 
have higher uh than the smooth wind but reduce the 
overall mass- loss rate. The case of clouds confined by 
a hot wind produces higher gas densities in the outflow 
than do unconfined clumps. Both models suggest that 
clumps may lower the inferred mass-loss rate by up to 
a factor of 10. We conclude that the mass- loss rates 
in warm outflows from 32016857 and TKRS 4389 are at 
least 30-50 M© yr~^and 35-40 M© yr~^, respectively. In 
this multi-phase model of the wind, we find the mass flux 
in the warm phase alone is comparable to the SFR. 

6. IMPLICATIONS 

One of the major goals of feedback studies is to mea- 
sure the rate at which galactic winds remove gas from 
star-forming regions. Our empirical understanding of the 
relative mass flux originated with the study of starburst 
galaxies nearby enough to directly detect their winds in 
X-ray emission. The mass flux in both the > 10^ K (Mar- 
tin et al. 2002; Strickland & Heckman 2009) and 10"^ K 
(Lehnert & Heckman 1996; Martin 1999) phases of the 
winds was shown to be of order the SFR. At higher red- 
shifts where winds presumably had an even larger impact 
on galaxy growth, the hot phase can rarely be directly ob- 
served; but unique observations of the warm-hot plasma 
now indicate it may contain 10 to 150 times more mass 
than the cold gas (Tripp et al. 2011). 

In this paper we demonstrated how the resonance emis- 
sion scattered by the warm gas in galactic winds can be 
used to measure the mass flux of low-ionization gas. Our 
discussion focused on the z = 0.9392 galaxy 32016857 
because we detected Mg II emission that is significantly 
more extended than the nebular [O II] emission and 
marginally more extended than the near-UV continuum 
emission. This halo emission almost certainly results 



from photon scattering because the Doppler shift of the 
Mg II emission is constant along the slit and does not 
share the prominent velocity gradient seen [O II] emis- 
sion. The halo gas appears to be an outflow based on 
the P-Cygni line profile of Mg II in the integrated spec- 
trum and the blueshift of the Fe II resonance absorption 
lines. This outflow must be collimated (i.e., an opening 
angle Q < 47r) because (1) the Mg II emission is only ex- 
tended in one direction (to the east) along the slit and (2) 
the optical depths of the Mg II emission and absorption 
components differ, indicating departure from spherical 
symmetry. 

Resolving the scattered emission improves our un- 
derstanding of the mass-loss rate in the low-ionization 
phase of the outflow. In the Sobolev approximation, the 
ionic density and the radial velocity gradient determine 
the scattering optical depth along a sightline through 
the outflow. Based on the Doppler shift of the low- 
ionization (Fe II) absorption and the spatial extent of 
the Mg II emission in the 32016857 spectrum, we esti- 
mate an average velocity gradient of 230 km s~^ over 
12 to 18 kpc for the outflow. Assuming we have de- 
tected the surface of unity scattering optical depth, an 
ansatz consistent with the doublet ratio of the Mg II 
emission, we directly calculate the density of Mg+ ions, 
n(Mg II) ^ (0.89 - 1.4) x 10"^ cm"^, in the outflow weh 
beyond the star-forming region. 

The implied mass flux depends on assumptions about 
metallicity, ionization fraction, depletion, and opening 
angle. We adopted a typical opening angle derived from 
the fraction of spectra of z ^ 1 star-forming galaxies 
showing blueshifted resonance absorption (Martin et al. 
2012) and a metallicity that ensures consistency with 
the mass - metallicity relation at similar redshift (Za- 
hid et al. 2011). The depletion of Mg onto grains is less 
constrained, but we argue that the hkely values span 
only a factor of four based on the well-measured proper- 
ties of Galactic clouds (Savage & Sembach 1996). The 
most uncertain parameter was the ionization correction, 
and we used photoionization modeling to show that the 
relation between the ionization fraction and the ioniza- 
tion parameter is nearly linear in the physical regime of 
interest. This calculation strongly suggests that most 
of the Mg in the outflow is Mg^ and therefore not di- 
rectly detected by our spectra. Following these phys- 
ically plausible arguments, the mass flux in a smooth 
outflow is approximately 330 - 500 Mq yr~^ in 32016857, 
where the estimated uncertainty is a factor of four and 
the SFR ^ 80 M© yr~^. We reached very similar conclu- 
sions for TKRS 4389 outflow using analogous arguments. 

Since we expect winds to be multi-phase, we further ar- 
gued that the dumpiness of the warm gas must be taken 
into account to accurately characterize the systematic 
uncertainty in the mass-loss rate. In an Appendix, we 
demonstrated that smaller clouds result in higher den- 
sities and lower mass-loss rates (relative to the smooth 
outflow). The requirement that the clouds survive out 
to at least the observed distance of the scattering halo, 
however, places a physical limit on the minimum size of 
the clouds; and we show that outflow models with clouds 
reduce the inferred mass-loss rate by a factor of ten rel- 
ative to the smooth wind model. Accounting for clouds 
lowers the inferred mass-loss rate in the warm outflow 
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but allows additional mass- loss in the hot phase. For 
32016857 and TKRS 4389, we reach a reasonably robust 
conclusion that the outflow rates in the warm phase are 
at least comparable to the SFRs of the galaxies and could 
be an order of magnitude higher in the absence of a hot 
phase. 

The mass loading factor, defined as r^ = M/SFR, 
serves a critical function in models of galaxy evolution. 
It scales inversely with the halo velocity dispersion, <j, 
for momentum-driven winds, but for energy-driven winds 
the relationship steepens to r^ oc a"^. This slope affects 
how quickly the faint-end of the galaxy mass function 
fiattens (Creasey et al. 2013; Puchwein & Springel 2013) 
as well as the slope of the mass-metallicity relation (Dave 
et al. 2012; Shen et al. 2012). To fit the galaxy mass func- 
tion and mass - metallicity relation, the mass flux asso- 
ciated with a given star formation rate needs to increase 
as galaxy mass decreases. The typical normalization in 
numerical simulations sets the mass-loss rate to roughly 
the SFR at some flducial mass. The relative mass frac- 
tion at widely different temperatures in these winds is 
not yet a robust theoretical prediction. Simulations fo- 
cusing on the winds rather than cosmological structure 
have just begun to include pressure from radiation and 
cosmic rays which analytic calculations suggest are im- 
portant (Murray et al. 2005, 2011; Breitschwerdt et al. 
1991; Everett et al. 2008; Socrates et al. 2008). At least in 
one recent simulation, Hopkins et al. (2012) find that the 
warm, 2000 < T{K) < 10^, phase dominates the mass 
transport in outfiows from SMC-like dwarf galaxies and 
remains comparable to the mass fiux in the hot wind in a 
model chosen to match the properties of a high-redshift 
starburst. 

At this time, TKRS 4389 (Rubin et al. 2011) and 
32016857 remain the only clear examples of Mg II scat- 
tered by outfiows. Among 145 LRIS spectra of 0.4 < 
z < lA galaxies, including 22 with very prominent Mg II 
emission, we found only three examples of spatially re- 
solved Mg II emission. Unlike 32016857, the line emission 
resolved around 32010773 and 12019973 is not necessar- 
ily scattered; it could come directly from H II regions 
as it is no more extended than the nebular [O II] emis- 
sion and not clearly distinct kinematically. We suggest, 
in spite of the large number of null detections, that dif- 
fuse Mg II emitting halos may be a common property of 
star- forming galaxies at intermediate redshifts. Approx- 
imately 15% of our sample have prominent Mg II emis- 
sion (Kornei et al. 2013). A larger fraction of the spectra 
show blueshifted resonance absorption lines indicative of 
outfiowing gas (Martin et al. 2012). The resonance emis- 
sion scattered by these outfiows could be attenuated by 
the LRIS slit or destroyed by dust absorption, a scenario 
supported by the observed tendency for Mg II emitters 
to have bluer than average U — B color (Martin et al. 
2012; Kornei et al. 2013). Based on the properties of 
TKRS 4389 and 32016857 - namely their high sSFRs 
and unusually high B-band luminosities for their U — B 
colors, we suggest that the scattered emission is brighter 
around these galaxies than it is near more typical Mg II 
emitters. 

While much of their halo emission may simply lie below 
current detection limits, intermediate redshift galaxies 
may prove an ideal environment for pinning down mass- 



loss rates. Outfiows were more common at z ~ 1 than 
today due to the higher cosmic SFR (Hopkins & Beacom 
2006) and smaller galaxy sizes (Simard et al. 1999; Lilly 
et al. 1998). Most of the star- forming galaxies have de- 
veloped well-defined disks by this time. The orientation 
of these galactic disks relative to halo sightlines probed 
by background quasars (or galaxies) provides a valuable 
means to distinguish intervening absorption created by 
outfiows from other sources such as infalHng gas (Bor- 
doloi et al. 2011). The relative metallicity of the galactic 
disk and halo clouds have also been used to identify ac- 
creting gas (Ribaudo et al. 2011). Scattered emission 
from an outfiow becomes easier to detect as the physi- 
cal area subtended by a fixed spectroscopic aperture in- 
creases, and the angular diameter distance increases with 
redshift out to z ^ 1.6. Since cosmic expansion reduces 
the observed surface brightness emission as (1 + ^)~^, 
the brightest scattered halo emission may be found at 
redshifts somewhat less than the maximum angular di- 
ameter distance. In the future, mapping the Mg II line 
emission scattered by circumgalactic gas should, in prin- 
ciple, provide a means to directly measure rj in galaxies 
over a broad mass range. The results from the first two 
detections of scattered Mg II emission strongly suggest 
very substantial mass-loss in the outfiows from interme- 
diate redshift galaxies. 
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TABLE 1 

Galaxy Properties 



Property 


32016857 


RA 


23:29:24.954 


DEC 


+00:07:05.85 


Redshift^ 


0.939196 


B (AB mag) 


22.49 


Mb - 5 log /i 70 (AB mag)^ 


-21.59 


U -B^ 


0.44 


logM*/M0^ 


9.82 


Vi ( km s-i)^ 


-91 ±15 


Vdop ( km s-i)^ 


-227 ±82 


\og[NDop{Fe+)Cf] ( cm-2)- 


14.41 - 15.53 


SFR ( M0 yr-i)d 


78 


sSFR (yr-i) 


1.2 X 10-8 



^From Martin et al. (2012). 

^From Willmer et al. (2006). 

^From Bundy et al. (2006). 
After scaling the flux of the red LRIS spec- 
trum to match the DEEP2 magnitude of Iab = 
21.70, we measured an integrated [O II] line flux 
of 6.25 X 10" ergs s" cm" , corresponding to 
an uncorrected (for extinction) line luminosity of 
2.79 X 10^^ erg s-^ To estimate the SFR, we 
use the average ratio of [O II] to Ha luminosity 
(0.40) measured by Moustakas et al. (2006) among 
nearby galaxies with Mb ~ —21.6. By adopting 
the mean ratio for galaxies with luminosities sim- 
ilar to 32016857, this calibration takes the strong 
correlation between line ratio and reddening into 
account. Finally, dividing the SFR obtained from 
the (Kennicutt 1998) Ho; calibration by 1.78, we 
obtain a SFR of 78 M© yr"^ for a Chabrier IMF 
from 0.1 to 100 Mq. 



TABLE 2 

Spectral Line Measurements 



Transition 




Ao 


Wr{X) 








(A) 


(A) 


Fell* 




2365.56 




-0.47 ±0.34 


Fell* 




2381.49 






Fell* 




2396.26 




-0.78 ±0.32 


Fell* 




2612.65 




0.31 ±0.33 


Fell* 




2626.45 




-1.41 ±0.32 


Fell* 




2632.77 




-0.32 ±0.24 


Fell 




2260.21 




<0.6 


Fell 




2344.21 




1.41 ±0.22 


Fell 




2374.46 




0.88 ±0.19 


Fell 




2382.77 




0.99 ±0.18 


Fell 




2586.65 




1.84 ±0.22 


Fell 




2600.17 




2.07 ±0.24 


Mgll 




2796.35 




0.34 ±0.39 


Mgll 




" 




-2.13 ±0.29^ 


Mgll 




2803.53 




0.85 ±0.22 


Mgll 




" 




-1.16 ±0.27^ 


CIIIJ 




1908.73 




-1.57 ±0.47 


CII] 




2326^ 




-2.06 ±0.37 


NeV] 




3426.98 




< 0.84 


Oil] 




3727.09, 


3729.88 


-128.2 ±3.1 


Nelll], 


He I 


3869.85, 


3868.63 


-14.16 ±0.95 


H8, He I 




3890.15, 


3889.75 


-6.83 ±0.63 


H I, [Ne 


III], He I 


3971.19, 


3968.58, 3965.85 


-8.73 ±1.03 


RS 




4102.90 




-9.55 ±0.70 


H7 




4341.69 




-16.68 ±1.36 



^Blend of A2324.22, A2325.41, A2327.65, A2328.84 

"In the integrated spectrum, the observed flux in Mg II emission is 
F(A2796) ^ 1.5 x 10"^^ ergs s"^ cm^^ and is F(A2803) ?^ 1.0 X 
10"^^ ergs s"^ cm"^ to roughly 20% accuracy. 
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APPENDIX 

A. MASS-LOSS RATE FOR A CLUMPY WIND 

The scattering optical depth of the A2796 Hne is 

T2796 = ^MgIlCr2796(Av) I. (Al) 

The cross-section of the A2796 transition, with units of cm^, is 

a2796(At;) = — /^, (A2) 

where Aq is the wavelength of the transition, e and rUe are the charge and mass of an electron, and / = 0.62 is the 
oscillator strength of the transition. The quantity Av is the velocity range over which the ions are distributed; for 
example, where the gas exhibits only thermal motions, Av = ^th, so that the velocity spread is simply the thermal 
velocity spread. The velocity spread may be much larger than ^'th if the gas exhibits bulk motions. 

In fact the expression (Al) for the optical depth depends on two unknown quantities, a characteristic length scale / 
associated with the outflow, and A'^. We scale Av to the ion thermal velocity '^th ^ 2(T/10^K)^/^ km s~^of a Mg^ 
ion, and introduce 

°V^ « 2.6 X 10-- f-^ V''%m^. (A3) 



^"- = ^^^-2.6xl0--(-^^ "^^ 

rrieC vth VIO J^. 



Using these, eqn. (Al) becomes 



T2796 = ^Mgll ^^2796^ ^ = ^Mgll ^^3795 ^eff, (A4) 



where we have introduced the effective scattering length 

W = ^l (A5) 

Av 

Specifying a physical model for the outflowing gas determines /eff- We will scale to a smooth wind model, in which 
the effective scattering length is the Sobolev length 

where dv/dr ^ v{r)/rs is the velocity gradient in the wind at the point where r2796 = 1. 

For our galaxy, taking the largest radius at impact parameter h where scattered emission is detected, 

A.l. Physical wind models 

We consider three models for the outflow, a smooth T ^ 10^ K wind, an outflow consisting of unconfined T ^ 10^ K 
clouds or clumps, and a two-phase hot {Th ~ 10^ K and smooth) together with a cold (T '^ 10^ K) and clumpy outflow, 
in which the cold clouds are confined by thermal pressure of the hot gas. The hot gas in the last model is naturally 
provided by supernovae exploding in the galaxy, at a rate proportional to the star formation rate (about one supernova 
per century for a star formation rate of one solar mass per year). 

A. 1.1. The smooth wind model 
In a smooth wind the characteristic velocity spread is Av = Vth and the characteristic length / is the Sobolev length 

W = feob = -TTT' (^^) 

dv/dr 



The logarithmic derivative of the wind velocity 



r(r).^^ (A9) 

V dr 



is of order unity in the wind acceleration region, so 

Zsob = r-1 (^) r « r. (AlO) 
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A. 1.2. Unconfined clouds 

Our second model assumes that the outflow consists of a large number of unconfined clouds. If the clouds are 
unconfined, e.g., if their gas pressure exceeds that of any hot outfiow, then their size along the direction of acceleration 
g is 

/ = ^^ (^) r^lOOpc, (All) 

where Cg ~ 10 km s~^ is the sound speed of the T ^ 10^ K gas in the clouds, and Vc is the circular velocity of the host 
galaxy, which is of order the outfiow velocity. The effective size of the cloud is 



eff 



mm(^- <*-' 



where we approximate Vc ^ v^ the outfiow velocity. 

The velocity spread A'^ ^ '^th if there are no bulk fiows in the clouds, in which case leff = Isoh- In the more likely 
case that there are transverse bulk fiows (expansions, in particular), then Av, could be of order c^. In that case, 

Q 

^eff ^ — feob ^ 0.1/sob- (A13) 

A. 1.3. Confined clouds 

Our third model assumes that the cold gas consists of clouds confined by hot gas. With a star formation rate of 
^ 8OM0yr~^, and an expected supernova rate of nearly one per year, it is likely that the galaxy possesses a hot 
{Th ^ 10^ K) gas component. Since the sound speed of this hot gas exceeds the circular velocity of the galaxy by a 
large amount, it will fiow out with a bulk velocity of a few times its sound speed, Ch ~ 1000 km s~^. This gas will act 
to pressurize the cold gas clouds, to accelerate them, and to disrupt them. 

A priori, the cold clouds could have any size Id and a range of internal velocities A'^;. However, we argue that there 
is a well defined lower limit to Id^ as follows. 

Both simulations (Klein et al. 1994) and experiments (Hansen et al. 2007) demonstrate that cold clouds embedded 
in such hot outfiows are destroyed on a few tens of cloud crushing times. The cloud crushing time is defined as (Klein 
et al. 1994) 

^) ^, (A14) 

rihj Vs 

where ric = n^f^ and Uh are the number densities of the cold and hot gas, and Vs ~ Vh is the velocity of the shock in 
the hot fiow. 
The clouds must live long enough to reach i?r=i ^ 11.4 kpc, which allows us to set a lower limit to Id'. 

30ree > {^^^f^\ , (A15) 

leading to a cloud length 

We have scaled the clouds survival time to 30 cloud crushing times. Scaling to a hot wind with a mass loss rate similar 
to the star formation rate, and assuming that A'^ ~ '^line, 

l^^ ^ 170 ( ^ ^ ( ^^ :.] PC. (A17) 

Vll-4kpcy V2 X lO-^cm-V ^ ^ ^ 

As with the unconfined cloud model, this length is similar to the Sobolev length associated with a smooth wind. 
The characteristic velocity spread Av ~ viine; it cannot be larger than the observed linewidth. Then 

W/=(^)^c^^l.Tpc. (A18) 

If the gas being removed from the cloud is over- ionized (so that all the Mgll ions are removed) before it is accelerated, 
Av could be as small as ^th, but as we will see this would lead to a large mass loss rate in the cold component. 

A. 2. Mass-loss rates 
For a smooth outfiow, the continuity equation 

M^ = nfcfhr'^nHV (A19) 

gives Eqn. 12. Note that r]Mg^ dMg^ and ey[g+ all enter to the one half power, so that the estimate for the mass loss rate 
is only moderately sensitive to these quantities. Observations indicate, however, that the warm outfiows are clumpy 
(Martin & Bouche 2009; Bordoloi et al. 2012). Taking /c ~ 1, where fc describes the fraction of the outfiow cone 
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covered by clouds as seen from the galaxy, places an upper limit to the mass loss rate in warm (10^ K) gas for a fixed 
density, ffinH- 
For illustration, we consider a bi-conical outflow of half-opening angle Ohc such that the outflow subtends a solid 

n 

angle VL = Vti^c = ^tt Jq ^"^ sin OdO. In light of the limited observational constraints on the orientation of the outflows in 
32016857 and TKRS 4389, the symmetry axis of the cone can be taken perpendicular to the observer's sightline to 
illustrate the impact of dumpiness on the density and mass-loss rate. 

A.2.1. Clumpy outflows 

The mass of a single cloud is 

47r 
rrici = -—llimnH,ci' (A20) 

Over the velocity width of the emission line, Avei, clouds at the projected radius where r2796(^) ^ 1 have to reflect a 
large fraction of the continuum. To find the number of clouds at this impact parameter which have projected velocities 
consistent with the linewidth, consider an annulus of width Ab on the sky. Over the velocity width of a single cloud, 
Avci, the clouds must cover a projected area 20b Ab. In addition, the clouds have to cover the range of velocity seen 
in the emission line. Together these two constraints give a lower limit to the number of clumps or clouds given by 

--/.■(^)(^)(i)(!?). 

where fy is the fraction of the linewidth covered by clouds. The total mass in clouds is Md = Ncirricu oi" 

"-Kv)4ts)(u)(!?)-f'>-— 



The mass loss rate is calculated using the crossing time, 

Ab Ab Ab 



r 



Vxib) '^sin^ V b 
We find 



2. ,^ ^ 2 flci\ fb\ fAvei\ 



(A23) 



M. = mr^v . 4^.,/.nH,. • - • ^f J ^- J ^^ J • (A24) 

Using equation (A19), we scale the mass- loss rate in the clumpy outfiow to the that of the smooth outfiow, 

--(l:)(f)^(^)-(l)(^)©Xt;)- 

It is convenient to scale the cloud density to a smooth wind: 

nH,ci = nH{Soh)(^-^) , (A26) 



where 

tQ K^ h^^Mgn(Sob)^ 
uh [Soh) = -f 

VMgCiMg^Mg+ 



(A27) 



For a clumpy outfiow, the effective length scale for absorption over the Hnewidth contributed by a single clump is 
simply related to the properties of a cloud, 

W/ ^ (^) Ici- (A28) 



\Av 
The density of a clump or cloud in the wind will exceed that of the smooth wind by a factor 

T X 1/2 / A X 1/2 

/Sob\ AVcl^ 



nH^ci = nH\^] ^^ . (A29) 

\lcl J \ Vth J 

A. 2. 2. Comparison of Clumpy and Smooth Outflows 
Using equations (A25) and (A29), 

The square-root dependence of the mass loss rate on the cloud size Id indicates this estimate is fairly robust. To 
illustrate how the clouds affect the mass-loss estimate, we take r = (1.0 — 1.5)6 ^ 6, 46>5c/^6c ^ 1, and fv/ fc ^ 1- 



For an unconfined clumpy outflow, with Id 



0.1 /sob and A'^ ^ c^, the mass loss rate 



Mci ^ 0.10 



280 



/ 10 y^' 



^Vth 



1/2 



hob 

90 



1/2 



1/2 



-^1 M. 
r 
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(A31) 



is about one tenth the smooth wind rate. The mass- loss is also lower when we assume hot gas conflnes the clouds. For 
a confined clumpy wind, with Id ^ Isoh but Av ^ vnne^ the mass loss rate is 



Md ^ 0.06 



Avd\ /280 Y^' 
280 y [AvdJ 



Avth 



1/2 



hob 

90 



1/2 



1/2 



12 



M, 



(A32) 



Allowing for very small cloud sizes Id^ as might be the case in a hot-gas confined clumpy outfiow, indicates that the 
mass loss rate can up to a factor of approximately 10 times lower than that estimated under the smooth outfiow 
assumption. 

Whether small clouds can survive the acceleration process remains open to some debate. The fiducial values of 9 and 
90 pc used here lie at the low end of the size ranges recently estimated from photoionization modeling of circumgalactic 
gas clouds in both nearby galaxies (Stocke et al. 2013) and intermediate redshift absorption-line systems (Meiring et al. 
2013). These calculations offer a physical argument nonetheless for the maximum amount by which the smooth outfiow 
model may overestimate the mass-loss rate. They suggest that the mass-loss rate in the warm outfiow could be just 
10% of the smooth wind mass-loss rate. If such clouds are present in the wind, then the hot wind will of course 
contribute additional mass fiux. 
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